Inelastic Electron Lifetime in Disordered Mesoscopic Systems 



Ya. M. Blanter* 

Institut fiir Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, 

Germany 
(February 6, 2008) 



Abstract 

The inelastic quasiparticle lifetime due to the electron-electron interaction 
(out-scattering time in the kinetic equation formalism) is calculated for finite 
metallic diffusive systems (quantum dots) in the whole range of parameters. 
Both cases of "continuous" (the inelastic level broadening much exceeds the 
mean level spacing) and "discrete" spectrum are analyzed. In particular, 
crossover between one- and zero-dimensional regimes is studied in detail. In 
the case of continuous spectrum the out-scattering time is shown to be the 
same as the inelastic time entering expressions for universal conductance fluc- 
tuations and persistent currents. It is also found to be shorter than the phase- 
breaking time in two- and one-dimensional systems, while in zero-dimensional 
systems these two times coincide. In the case of discrete spectrum for small 
enough systems a universal behavior of the scattering time is obtained. For 
temperatures below the mean level spacing the out-scattering rate is shown 
to be vanishingly small. 
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Typeset using REVT^ 



I. INTRODUCTION 



As is well known, inelastic electron scattering plays an important role in various phe- 
nomena in disordered metallic systems (see the extensive review by Altshuler and Aronov 
HI]). It is enough to mention that it is responsible for the weak localization correction to 
the conductivity. It proved also to be important in mesoscopic phenomena in these sys- 
tems. In particular, inelastic scattering governs the temperature dependence and crossovers 
between different dimensionalities for universal conductance fluctuations (UCF's) and 
similar problems, such as persistent currents (see e.g. 0]) and correlators of persistent cur- 
rents within the grand canonical ensemble (GCE) (see e.g. [Q). Recent estimations ^ have 
shown that the conductance of mesoscopic systems in the case of a discrete spectrum, at least 
within the canonical ensemble (CE), is also governed by the inelastic electron scattering. 

The electron-electron interaction enters these problems through the inelastic scattering 
time. Extensive investigations carried out about 10 years ago showed that at least two 
relevant electron-electron scattering times exist : (i) the out-scattering time Tout appearing 
in the kinetic equation formalism ; it has the meaning of an inverse frequency of inelastic 
collisions; (ii) the phase-breaking time (see also [0]), responsible, in particular, for 

the weak localization correction and for the quasi-particle decay within the Fermi liquid 
theory. These two times has been shown to coincide in three-dimensional systems {t^ ~ 
Tout oc T^^/^). The phase-breaking time is very well studied also in infinite 2D (r^ oc T^^) 
and quasi- ID (r^ oc T~^/^) systems [0]. However, the out-scattering rate in low-dimensional 
systems 7(e, T) = ToJt have caused some controversy. It was studied for zero temperature by 
Altshuler and Aronov (see [0]) who obtained the result 7(e, T = 0) oc e'^/^, d = 1 (QID) or 
d = 2 (2D). This implies, in particular, that the Fermi liquid theory is violated in quasi-lD 
systems close enough to the Fermi level. An attempt to include the finite temperature makes 
the situation even worse: the out-scattering rate diverges for low enough energies e <^T . In 
two-dimensional systems attempts to cure this singularity have been made [jl2|-0, leading 
to different results. 

Below this problem is revisited. The out-scattering rate 7(e, T) is considered for T ^ e. 
In principle, one should study the two-parameter problem, but in order not to make the 
expressions too cumbersome and not to consider a huge number of parameter ranges, e is 
just put to be equal to zero [0. It is shown that the out-scattering time definitely does 
not coincide with the phase-breaking one. For the 2D case (Sec. H) just the result by Refs. 



l^Hl4| 7 oc TlnT is recovered. In the quasi-lD system (Sec. |TD new results are obtained 



for large temperatures the out-scattering time shows the same temperature dependence as 
the phase-breaking one (7 oc T^/^) but is parametrically shorter than the latter; moreover, 
for low enough temperatures (but still in the metallic regime) the out-scattering rate exceeds 
the temperature, and exhibits another temperature dependence 7 oc T'^^^. 

Now we return to the effect of inelastic scattering effect on the mesoscopic phenomena in 
disordered systems, in particular UCF's. It is rather clear (and it is shown once more below) 
that the inelastic scattering time entering the UCF problem is essentially the out-scattering 
time. (Note that because of the controversy mentioned above it it is believed also to coincide 
with the phase-breaking time, that is definitely not the case for 2D and QID systems). So 
one needs to investigate the out-scattering rate in finite systems. The problem of electron 
lifetimes in finite systems has not been addressed so far, except for the only paper by Sivan, 
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Imry and Aronov [|T8| (to be discussed below, see Sec. ^) . 

The introduction of the finite size of the system leads to the appearance of two charac- 
teristic energies: the Thouless energy E^. = D / L'^ (here D and L are the diffusion coefficient 
and the largest size of the system, respectively) and the mean level spacing A. The small 
parameter of the standard perturbation theory (in particular, the diagram technique) is 
7/A, and respectively two cases should be distinguished: "continuous spectrum", 7 ^ A, 
and "discrete spectrum", 7 ^ A. 

In the three-dimensional systems (e.g. cube of the size L) the situation is rather sim- 
ple: for T ^ Ec the inelastic level broadening 7 much exceeds the mean level spacing 7 
(continuous spectrum), and the results actually do not depend on L. However, for T <^ Ec 
one obtains 7 ^ A, and the crossover to the zero- dimensional situation occurs. So in the 
case of 3D - OD crossover the distinction between 3D and OD situations is the same as 
one between continuous and discrete spectra (for the exact formulation of this point see 
Sec. fVp. However, it is not the case for 2D - OD and QID - OD crossovers. The out- 
scattering and the phase-breaking rates are calculated in Sees. || (2D - OD crossover) and 
pn| (QID - OD crossover) in the case of the continuous spectrum by means of the diagram 
technique. It is shown, in particular, that the spectrum for T ^ Ec is always continuous, 
and for zero-dimensional situation (the result depends on L) the two inelastic times coincide: 

Tout ~ T^- 

In the case of discrete spectrum 7 <^ A the perturbation theory is no longer valid, and 
the random matrix theory (RMT) should be involved. Two cases of "small" L <^ ppl'^ 
and "large" L 3> ppl"^ dots should be distinguished. The temperature dependence of the 
out-scattering time is obtained 7 ~ T'^A/E^ for A < T < in "large" dots, while "small" 
dots exhibit an universal behavior 7 ~ T'^/ep. For the temperatures below the mean level 
spacing T ^ A in all cases the scattering rate 7 is found to be vanishingly small (Sec. |rV|). 

Below diffusive metallic systems are considered. "Diffusive" implies the following relation 
between the characteristic lengths of the problem: / ^ Li with / and Lj being the elastic 
mean free path and linear sizes of the system respectively. For definiteness a rectangular 
sample with sizes = L > Ly > = a is chosen; h = 1. "Metallic" means that (i) 
the disorder is weak, ppl ^ 1, pp being the Fermi momentum, and (ii) in quasi-lD system 
the localization length ^ ~ l{ppaY is larger than either L or the phase-breaking length 
{Dt^Y/"^. The closed systems are of primary interest. This means, that the coupling to the 
environment is assumed to be weak, namely -C A, Tc being a characteristic time of the 
electron escape through the attached leads. For open systems the single-particle states are 
smeared; the magnitude of this smearing is of order 7^ ~ min{i?c, t~^}, and respectively 
one is interested in the energy (temperature) range e ^ 7c. The results obtained below are 
valid, however, for the open systems also provided 7 ^ 7c (see the discussion in Sec. 

II. CONTINUOUS SPECTRUM: GENERAL EXPRESSIONS AND 2D CASE 

A. Generalities 



In the case of continuous spectrum 7 ^ A the standard impurity diagram technique |T9 



may be involved. We use below the diffusion-cooperon approximation For T ^ A one 
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can not expect any difference between GCE and CE; as it will be shown below the spectrum 
is always continuous for T ^ ^ A, so all calculations are performed for the more simple 
GCE case. 

It is well known that in 3D system the phase-breaking time 

-1 ^ rT/n^3/2 -1 



{T/DY^W (1) 



coincides with the out-scattering time and is governed by large momentum transfer Dq^ ~ T. 
Here D = t^/rd is the diffusion coefficient {d is the dimensionality, presently d = 3) while 

= itipf/tt^ and r are the density of states and the elastic scattering time respectively. 
So 3D case is not interesting for us. Crossover to other dimensionalities occurs when one 
of the sizes is under (Z^/T)^/^. Below low-dimensional systems are studied; the inequality 
T D/a? \s assumed to be satisfied. 

In the UCF theory and related problems the inelastic scattering time appears as a result 
of calculation of the "UCF diffusion propagator" [Q]. It differs from the "true" diffusion 
propagator |l|] by the absence of electron-electron interaction lines connecting two electron 
propagators. A "true" density-density correlation function is not affected by the electron- 
electron interaction due to the Ward identity, that is responsible for the cancelation of 
interaction lines from the vertex corrections and from the corrections to electron propagators. 
However, if some diagrams are omitted, the Ward identity does not take place any more, 
and hence the electron-electron interaction renormalizes the UCF diffusion propagator. In 
the absence of magnetic field the "UCF cooperon" is given by exactly the same expressions 
as the UCF diffusion propagator. 

The diagram equation in the lowest order of the electron-electron interaction for the 
diffusion propagator is shown on Fig. 1 []8|,p!3[|. The shaded rectangle is the full diffusion 



propagator r(q, a;), the double dashed line is the bare diffusion propagator 

TXUdr^[\uj\ + Dq^) 

while the crossed rectangle is the block concerned with the electron-electron interaction 
J(q, cj), displayed for the UCF case on Fig. 2. The density of states in quasi-ii dimensional 
system Vd (to be referred below as v) is equal to Vd = i'3a^~'^; for pure 2D system i>2 = m/ir. 
The Matsubara diagram technique is used; it is assumed for definiteness that u + e > 0, 
e < 0, and after the analytical continuation iu ^ u, ie ^ e the propagator corresponds to 
the averaged product of retarded and advanced Green's functions {G^iu + e)G'^(e)). The 
equation can be readily solved to obtain: 

= 2i\ I ! n 2 I — V' (^) 

T[VT^[\uj\ + Dq-^ + 7) 



7 = -(ttz/t^) ^/(q,cj) 

So the inelastic rate 7 is given by the interaction block and is essentially the same as the 
out-scattering rate in the kinetic equation approach 0. 

Now we calculate the interaction block. In the arbitrary dimensionality one obtains for 
the diagrams of Fig. 2a, 2b, 2c and 2d respectively (note that for Fig.2a and Fig. 2b the 
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leading terms cancel and one has to perform an expansion over the small momenta and low 
frequencies) : 

/.--^^TEE^^^(-^ + ^^^^ (4) 



Here U (q, 1)) is a screened Coulomb interaction. 

B. 2D case: Infinite system 

As the Coulomb interaction has different forms in different dimensionalities , the ana- 
lytical continuation should be performed separately for d = 2 and d = 1. We start from the 
2D case with 



27re2 




\+Dq^ 
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fi| + Dkq 



[/(q,fi) = — (5) 



Here k = Aire'^i'^a for quasi-2D case and n = 27re^z/2 for 2D case. 

After the analytical continuation iuj —>■ u, ie ^ e one obtains for = 0, q = 0: 

g2 /" / Z Z c\ Z 

^=nL ^^^g^^^y dz [coth - - tanh ^ ^^^^^^^^ ^ ^^^^^^^ (6) 

It is important that the terms with the hyperbolic cotangent come from the quantities Ic 
and Id- 
It is easy to see that for non-zero temperature the integral diverges at 2; = 0. This 
singularity can be cured by the substitution of self-consistent expression for the diffusion 
propagator (H) to the formulae (|) instead of T^^\ The diffusion propagators entering the 
diagrams of Fig. 2a and 2b are "true" and are not renormalized by the electron-electron 
interaction, while those of diagrams Fig. 2c and 2d are renormalized. For T 3> e the 
main contribution to the integral comes from the region \z\ < 2T, and one obtains the 
self-consistent equation for 7 (that is supposed to be momentum- independent for a while): 

4e'^T f°° / \ /■^'^ dz 
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A straightforward calculation leads to an equation 




(8) 



where the last equality is the leading term of the solution. Substituting the explicit expres- 
sions for K, one finally obtains: 

^ / (T/2e^r) In {SDK^epTy/T) 2D 
' \ (37rT/4p|a/)ln(DK2(4p^a//37r)2T-i) Quasi-2D 



Note that this rate is higher (by a logarithmic factor) than the 2D phase-breaking rate |T0 



^ \ {pFay^\n{pj,la) Quasi-2D ^ ^ 



The result (^) is valid under the condition T 3> 7 only. In principle, this condition (unlike 
the condition of the Fermi-liquid theory validity Tr^ ^ 1) is violated in the exponentially 
small range of low temperatures. We consider this problem as purely academic, however, 
and do not write down expressions for the corresponding range. 

The expression (^ was derived first in Ref. [0 as an imaginary part of the impurity- 
averaged self-energy. Later this result was confirmed in Refs. 0,0] (Fukuyama and Abra- 



hams |TB[ have used the technique identical to that of the present paper). However it was 
believed to be the expression for the phase-breaking time in a 2D system. The origin of this 
confusion lies in the statement |13[ that the diagrams with the interaction between differ- 
ent electron lines (i.e. those not contributing to 7, but improtant for r^, in-scattering terms 
in the kinetic equation approach [§]) are small. However, the original calculation of the 
phase-breaking time [10|| allows one to separate these two contributions (with and without 



interaction lines between different Green's functions) explicitly, and the result is that these 
contributions are of the same order. So the difference between 7 and t7^ should be looked for 



in the diagrams omitted in Ref. . Hence the attempt to cure this discrepancy by the 
introduction of a small-momenta cutoff does not seem well justified. Indeed, the momenta 
range proposed to be cut off is exactly the region of relevant momenta in the integral, and 
the introduction of finite momenta does not change the result. So we argue that the result 
by Fukuyama and Abrahams, (^), is true; however, it describes not the phase-breaking time, 
but the out-scattering time (or the UCF inelastic time). We have shown that the inelastic 
time entering the UCF problem is essentially the out-scattering time and differs from the 
phase-breaking one. It will be seen in Sec. |TTl| that the difference is even more pronounced 
for the quasi-one- dimensional case. 



C. 2D case: Finite system 

Now we turn to the description of 2D - OD crossover. In the derivation of Eq. (|^) it 
was supposed implicitly that the system is infinite. Let us now consider 2D finite system. 
Under the conditions of diffusive regime L I, and in the case T ^ the equation (^ 
is still valid, however the integral should be understood as a sum over momenta with the 
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boundary conditions taken into account. So q = {tt/ L){nx,ny) with integers Ux and Uy 
allowed or not to be equal to zero, subject to the boundary conditions. It is important that 
in any case the q = mode does not contribute to the integral: for open systems (attached 
to the metallic junctions, i.e. in the UCF problem) this mode is forbidden by the boundary 
conditions, while for closed ones (e.g. in the problem of persistent currents) the conditions 
of the charge neutrality imply U{0,u!) = 0. So the difference between these two cases is not 
quite important, and the boundary conditions will be taken into account just by the cutting 
off the integral for small momenta q < n/L. 

The analysis of the equation (|^) then shows that for 7 ^ (EcDk'^Y^'^ this cutoff does 
not play any role and the result (^ is still valid. Note that this region corresponds to the 
2D case from the diffusion point of view. For the opposite case, that is realized under the 
conditions 

r(//L)epr, 2D 
/ ^ ^ \ {l/L){plal), Quasi-2D 

one obtains 

e^T , 2T 
In — — 



7=^ln^, (11) 



It is seen that this result also differs from Eq. (|]) by a logarithmic factor. 



III. QUASI-ID CASE AND QID OD CROSSOVER 

Let us now turn to the quasi- ID case. The screened Coulomb interaction is 

U{q, O) = e In \n\ + Dq^ + e^UiDq^ \n{q^a^)-^ ^^^^ 

It is convenient to denote the denominator of this expression as \Q\ + Dq^f{q) with /(g) = 
1 + Cln(g^a^)~\ C = e^vi. If, as usual, ~ f, the constant C occurs to be of order 
C ~ {ppaf > 1. 

Prior to the calculation of 7 for the quasi-lD system it is necessary to analyze the 
expression for the phase-breaking time [|1^] that can be conveniently rewritten as 



This result was obtained under some assumptions to be analyzed here. First, the Fermi 
liquid theory is valid only in the case Tr^ ^ 1. Second, the system is one-dimensional, i.e. 
Ec ~ D I L"^ <^ T^^ <^ T -C D/a?. Finally, it is supposed to be metallic, i.e the correlation 
length i = l{pFCif > min{L, {Dr^f^}. 

The last condition is violated in the localized regime (the shaded range on Fig. 3): 

Tt < (pFa)"^ l/L < (pFa)"^ 

and the first three are consequently summarized as (Fig. 3, ranges I, II, III, and IV): 
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{l/Lf{pFaf <^Tt<^ {l/af < 1 (14) 

The question is of much interest what happens if the condition Tr ^ {1/ L)^{ppaY (for 
1/ {Pfo) "^) is violated. The expression for the phase-breaking time in an infinite system 
(p!3|) is no more vahd there. However, it is possible to derive the expressions for this regime. 
This is done in the Appendix A. The result is 

i ~ (16) 

{l/Lf ^Tr^ {l/Lf{pFaf (16) 

This expression is valid until ^ T (range V), as is seen from Eq. (|TBp. As will be shown, 
in the whole range (0) the inequality 7 ^ A is satisfied, A = {v^V)'^ being the mean 
level spacing. So the spectrum is continuous and this range is subject to the perturbative 
analysis. Thus, one has a "true" quasi-lD region (p!^), transition QID - OD region (|I6|) 
and "true" OD region T <^ where the spectrum is discrete. In this Section both QID 
and transitional ranges are considered. It will be shown that a new splitting of the range in 
respect to the out-scattering time appears. 

Now it is possible to start the calculation of the out-scattering time. Consideration of the 
diagrams Fig. 2 leads to the divergent expression, and this divergence can be cured exactly 
in the same way as in 2D case. Extracting the main term for T ^ e, one obtains: 

Here a low-momentum cutoff is introduced in the momentum integral (see the discussion 
for 2D case). Equation ( [T7| ) contains all information about the scattering rate and should 
be solved in different limiting cases (Fig. 3). 

A (Range I). Let us assume D{t[ / L)"^ fiji / L) ^ 7 ^ T. In this case the cutoff is not 
important, and the essential momenta in the g-integration are Dq^f^q) < 7. The equation 
simplifies to a form: 

7 = ATeyCDqon, Dqlf{q^) = 7 
(the inequality C ^ 1 has been used). This yields the result: 
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TT 



(CD)-^/=^lni/3 



2/3' 



a2 



\4Te^J 



(18) 



This cumbersome expression can be simplified if one assumes again ~ 1;. Then, as the 
cube root of the logarithm is always a quantity of order unity, one obtains: 

^^T2/3r-i/3(p^a)-2/3 (19) 

It is seen that the initial assumptions D{7i/L)^f{'K/L) ^ 7 and 7 ^ T are satisfied only in 
the temperature range 

max{(pi.a)-2, {l/L)%pFa)^} <^Tt <^ {l/af 
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(Range I). In particular, this result is valid for the infinite quasi-lD system for temperatures 
above r~^(pi?a)~^. It should be emphasized also that in this case r^'y = (pfI)^^^ ^ 1, and 
consequently Tout = ^ t^- So the out-scattering time in this range of parameters shows 
the same temperature dependence as the phase-breaking time, but is much shorter than 
the latter. Nevertheless, the transition to the OD behavior in the problem of out-scattering 
times occurs before the corresponding transition in the dephasing problem. 

B (Ranges II, III). If one assumes D{ti / LY fiji / L) ^ 7, 7 ^ T, the straightforward calcu- 
lation leads to the equation as follows: 

A new relevant length scale appears, g, defined by Dqf{q) = 2T. Then two limiting cases 
should be distinguished: 

1. (Range II). q ^ L^^. The cut-off is again unimportant, and one obtains: 

2eT3/4 , CD 

7 



7rV2(CD)i/4 



In-^/^ (^) ~ r3/V-V4(p^a)-V2 (21) 



The second identity is the estimation, based upon assumption ~ f . All assumptions 
made for the derivation of Eq. ( pID occur to be consistent in the range II: 



max{ (pi^a) ^, (//L)^(pira)^} <^Tt <^ {ppct) ^ 

This result does not contain the length of the sample again and consequently describes 
the infinite quasi-lD systems for the temperatures just above the metal-insulator transition: 
{pFa)~'^ -C Tt <^ (pfo)''^. We have obtained a new temperature dependence, 7 oc T^^^, 
and the out-scattering time is now not only parametrically shorter than the phase-breaking 
one, but even 7 ^ T. As 7 has a meaning of a frequency of electron-electron collisions, the 
violation of the condition 7 <^ T has nothing to do with the violation of Fermi liquid theory 
(the latter is subject to another condition, Tr^ ^ 1). 

2. (Range III), q <^ L~^. Now the scattering rate introduced to the rhs of Eq. (^Of ) for 
self-consistency, is inessential, and one obtains 

This result is valid under conditions (range III): 

l/L < (pFa)~^ {ppa)'^ ^Tt<^ {l/L)\pFa)-^ 

Now the result contains the length of the sample, but still one has 7 ^ T and 7r^ ^ 1. 

C (Ranges IV + V). In the case D{'k / LY f{'K / L) 3> 7 the important region in the integral 
over 2; is 2; < Dq^ + 7, and again straightforwardly all dependence on the inelastic time 
cancels out in the rhs. The result is 



9 



[l/Lf <^Tr<^ {l/Lf{pFaY 

To summarize, we have described the crossover between quasi-one- dimensional and zero- 
dimensional behavior for the case of continuous spectrum. We have discovered 5 different 
parameter ranges (Fig. 3), that can be divided into three groups: 

1. I and II. Here both the dephasing and the out-scattering are purely quasi-one- 

dimensional, the results do not contain L. In the whole range the relation Tout 't<p 
holds. 

2. Ill and IV. The transitional region: the dephasing is still quasi-one-dimensional, while 

the out-scattering time depends on L and is of zero-dimensional nature. Still one has 

'''out ^ T<f>- 

3. The range V is truly zero dimensional, the expressions for the phase-breaking rate ([ISD 

and for the out-scattering rate (^31) coincide. However the spectrum still is continuous, 
i.e. r^^ ~ T~ut > A. 

The spectrum becomes discrete for the temperatures below the Thouless energy: T ^ E^, 
as it takes place in 3D case. 



IV. DISCRETE SPECTRUM 

The diagram technique used in Sees. |TT| is valid only in the case of continuous 
spectrum 7 ^ A. As was shown above, this condition in all cases is equivalent to another 
one T 3> Ec- So this perturbative approach is not suitable for the case of low temperatures 
(or, alternatively, small dots) T <^ E^. Fortunately, this is exactly the range where the 
random matrix theory (RMT) works quite well for the description of small disordered 



systems |21]. Note that usually RMT in the range of its applicability is equivalent to the zero- 



dimensional non-linear sigma model p2|. However, in the problem under consideration one 
should use the four-point correlation function which up to now was not derived by means of 
the supersymmetry method. So RMT (in spirit of the paper by Gor'kov and Eliashberg pip 
seems presently to be the only method for investigation of the electron-electron interaction 
in the non-perturbative regime. To avoid misunderstanding, we stress that the electron- 
electron interaction is taken into account perturbatively; however, the parameter A/7 is no 
more small and should be treated non-perturbatively. 

In order not to overburden our expressions, we consider in this Section the cubic sample 
with size L. The system is assumed to be charge neutral, and the results rely heavily on 
this fact. In principle, the methods used allow one to consider the samples of arbitrary 
geometry, and all expressions below containing the parameters A and E^. instead of L are 
valid for the general case (for the metallic regime). Also, we do not distinguish two types of 
the scattering rates from each other. To understand the results better, we extend the ranges 
of parameter to the clean case also. The exact result for the scattering rate in bulk 3D case 
reads as IHl 
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3l/2 ^ 3/2 2 

^^^'^ = 0) = ^r^^ + « - (24) 

The first term is exactly Eq. (p. It exceeds the second one for small enough temperatures, 
Tt <^ and in the previous sections it was implicitly supposed that this condition 

is satisfied. However, for Tt ^ (pf^)'^ the second term in the most important. This is just 
the result in the "clean" case: the mean free path does not enter the expressions. So we 
consider in this section the case L ^ / for arbitrary temperatures: both diffusive and clean 
limits. The result (|^) is obtained by the diagram technique and is valid consequently for 
7 ^ A only; in the diffusive limit this condition is equivalent to T ^ E^, while in the clean 
limit for T = Ec one has still 7 ^ A. The methods used below are valid for T <^ Ec, and so 
in the clean regime one has some overlap between the results to be obtained and Eq. (p4D. 
Below the out-scattering time 7"^ for both the GCE and CE cases is calculated. 

Prior to the the calculation it is necessary to analyze the paper by Sivan, Imry, and 



Aronov |T8|], devoted to the inelastic scattering rate in the finite systems. The authors of 
Ref. |]TB[ calculated exactly the same quantity as we do — the scattering rate as a function of 
energy 7(e). The methods they used are valid in the case of "continuous" spectrum 7 3> A 
where they have just reproduced the results (0). However, the formal application of these 
methods to the opposite case of "discrete" spectrum (7 ^ A) leads to some non-trivial 



results, that appear to be not well justified in the approach of Ref. [18|. Below we present 



a rigorous calculation, valid in the case 7 <C A, and the results in the range T ^ A (where 



level correlations are unimportant) exactly reproduce those of Ref. [jT8|. The reasons for this 



equivalence are still unclear. Also the effect of level correlations, which can not be taken 



into account by the methods of Ref. [|T8|, is authomatically incorporated in the approach 
developed below. 

A general expression for the scattering rate, based on the perturbation theory of the 
Coulomb interaction, is 

7A = 27r |(A,A2|f/(ri,r2)|Ai,A3) p5(eA + eA,-eA,-eA3)/A.(l-/Aj(l-/A3) (25) 

-^11-^2,-^3 



(A,A2|f/|Ai,A3) = J rfridr2f/(ri,r2)^I(ri)^A,(ri)^I,(r2)^A3(r2) 

Here |Ai) = ipXii^r) are exact single-particle states - single-electron states in the unique 
disorder realization in a quantum dot (note that all four states need to be different) — fx^ 
are Fermi distribution functions (the CE case needs a more delicate treatment though; see 
below) and U stands for the electron-electron interaction. The expression (|25|) should be 
averaged over the disorder realizations. 



A. Coulomb interaction in the restricted geometry 

The first thing is to calculate the screened Coulomb interaction in the restricted geometry. 
For the temperatures (energies) below the Thouless energy one can study the static screening 
only. The latter is a solution to the equation as follows: 

V2t/(ri, r2) = K''9iri)eiv2)U + Aire'^Siri - r2) (26) 
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K = (47re^z/)^/^ being the inverse Debye length, z/ = z/3. The functions 9{r) are equal to 
unity and zero inside and outside the sample respectively. This equation can not be solved 
exactly for a rectangular sample, and one should introduce some approximations. 

In order to make these approximations clear, we consider first Eq. ( PB]) in a simplest 
geometry, where it can be solved. Namely, if the sample occupies a half-space x > 0, the 
exact solution to Eq. (^) in the region a; > 0, x' > has the form: 

t/(ri,r2) = J ^^^^ exp [iqy{y - y') + iq^{z - z')] fg{x, x'), (27) 

fq{x, X') = — exp {—p\x — X'\) P___^ _|_ ^/-jj ^ p = (^g2 _|_ ^2-jl/2 j^28) 

Here the first term is essentially the screened Coulomb interaction in the continuous media; 
it is translationally invariant and for small screening length is proportional to 6{x — x'). 
The second term is due to the restricted geometry effects; it is important that it be essentially 
nonzero when both x and x' lie in a narrow layer along the boundary of the sample; the 
thickness of this layer is of order k~^. 

Now consider another problem. Let us impose the boundary condition in Eq. (p6|): 
U |x=o= (and respectively U is nonzero only when both x and x' lie inside the sample). 
Then the both 6'-factors are identically equal to unity, and we obtain: 

fn{x, x') = exp {—p\x — x'\) -\ exp [—p {x + x')] (29) 

p p 

If another boundary condition is imposed: dU/dx \x=o= 0, a solution is 

. r. 27re2 , , 

fg{x,x) = exp{~p\x-x\) (30) 

p 

It is seen that the results the problems ( pOf ) and (|30|) yield inside the sample (x > 0, x' > 
0) differ from exact ones (|28|) by the contribution that is nonzero only near the boundary 
of the sample. Now one should recollect that the initial problem requires only the matrix 
elements of the screened Coulomb potential, and so one has to consider the region inside 
the sample only, and the contribution of the boundary term is small in comparison with the 
main one by a factor {kL)~^ <C 1. So Eq. (^) may be replaced by a more simple one, with 
6'-factors set to be equal to unity and the boundary condition imposed. Moreover, since the 
range of the interaction is k~^, i.e. is extremely small, this result is not sensitive to either 
boundary conditions or the sample's geometry, and one may perform this replacement for 
our rectangular sample as well. 

So, turning to the case of the rectangular sample, one easily obtains: 

C/(r.r')^Wi:;^iW^ (31) 

Here (Pq{r) are the eigenfunctions of the Laplace operator with appropriate boundary con- 
ditions; the eigenvalues are equal to — g^. In the case of the specified cubic geometry one 
gets 

q= {n/L){nx,ny,n^), ni = 0,1,2,... 
If the system is charge neutral, the q = mode should be dropped in the summation. 



12 



B. Calculation of the matrix element 



Now we return to the equation (P^D- Actually the squared absolute value of the Coulomb 
interaction matrix element contains the product of eight single-particle eigenfunctions, and 
our statement is that this combination only weakly depends on the energies of these states, 
and therefore can be impurity averaged separately from others, energy dependent, factors, 
yielding the constant f/^. The results obtained below follow in fact from the quasi-classical 
approximation pT| , |23[ . 



One needs to average 

J(ri,r2,r3,r4) = (V^I(ri)^A,(ri)^A(r3)C(r3)V^A,(r2)V^A3(r2)V^A.(r4)^l3(r4)) (32) 
so that 

= /*,*.A<;r,C/(r,.r.)C/(r3.r,)J(r,.r„r3.r,) (33) 

In the particular case (pi]) f/(r, r') = v^^5{y — r') one obtains 

= j dvidv^Jivi, ri, rs, rg) (34) 

The quantity J contains the constant part, that corresponds to the calculation in the 
Gaussian ensemble, and coordinate-dependent contributions of higher modes. 

Gaussian ensemble. The constant part can be easily calculated. It is convenient 
to introduce the discrete notations. If one considers N electrons (A^ ~ ^f/^ ^ 1), and 
consequently splits the system over elementary volumes with positions Tj, then the values 
of each two eigenfunctions in each two elementary volumes are in the Gaussian ensemble 
independent except for the constraints due to the orthogonality and normalization conditions 
for these eigenfunctions. In the leading order in A^~^ one has: 



JG(ri,ri,r3,r3) 



V-\ ri = rg 

2N-^V~^, ri ^ rs 



2T/-4 ^ / ^ (35) 



V being the volume of the system. 

Hence the contribution to f/^ from the Gaussian ensemble is 

f/2 = y-'^V-'^N-'^ = A'^/N = aA^/ep (36) 

Here a is a numerical coefficient. It can be adjusted from the comparison with the clean 
limit in the overlap parameter range (see below); this gives a = 1/8. Note that the main 
contribution to the integral (|33D comes from the range where all four coordinates coincide; 
in the continuous notation, this means that the distance between these points is of order of 
the screening length. 

Higher modes. The contribution of the higher modes is concerned with the diffusion 
processes. In particular, the coordinate-dependent part of the equation ( p2D describes the 
diffusion of an electron from point ri to point r^, and the diffusion of another electron from 
point r2 to point r4. It is reasonable to assume that if point ri is far enough from point r2 
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(in the discrete terms, these two points he in different elementary volumes, or, alternatively, 
the distance between these points exceeds several interatomic distances), and point is far 
enough from point r^, this diffusion processes are independent: 

JHAf(ri,r2,r3,r4) = (V^^(ri)^Ai(ri)^A(r3)V^I,(r3))(^l2(r2)^A3(r2)^A2(r4)V'l3(r4)) (37) 

If, however, these pairs are close to each other (in particular, this is the case for the short- 
ranged interaction), additional contributions from the diffusion process ri —>■ r^, r2 —>■ 
appear, and the expression ( pTj) acquires a coefficient 2. 

The average of four eigenf unctions of the type (|37D can be calculated up to the terms of 
order g^"^, with g = Ec/A being the conductance. The result reads [pl| , p^ : 



{rx{r)^x'{r')Mr')rx'{r')) = iN^V) (W^r(r',t) - V-') dt (38) 

Here Wrir', t) is the probability to find an electron at the time moment t in the point r' if it 
was in the point r in the time moment t = 0. This probability obeys the diffusion equation: 

dW/dt = DVl,W, W{r, t = 0) = 6{r' - r) (39) 

Integrating Eq. (PU]) over the time variable and taking into account that for t — >■ oo the 
distribution tends to be uniform one: W = V~^, one obtains: 



L 



Wdt = + ^(Dg2)-V.(r)¥^,(r') (40) 



Finally, combining Eqs. (0), (0), (0), (pH]), and (^), one obtains for the contribution of 
the higher modes to the interaction matrix element: 



UIm = ^ (41) 



The constant S is given by 



If one rewrites the contribution of the higher modes as U'^ ~ A'^g^'^, g = Ec/A being 
the conductance, it is easily seen that it has the same contribution as that of the Gaussian 
ensemble, except for the factor A^~^ being replaced with g~'^. Consequently one obtains 
^g/Uhm ~ q'^/N 123 ■ large sample, L ^ Pf^"^, "we obtain g"^ -C A^, and consequently 



the contribution of higher modes is the leading one. In the opposite case, L -C Pf^'^, the 
Gaussian ensemble produces the leading term. The latter is universal, i.e. does not contain 
any information about the disorder. It is shown below that this contribution leads to the 
clean-limit result 7 T^/e^. 

To summarize, we have obtained the following expression for the averaged matrix element 
of the Coulomb interaction: 



U'=\ (A,A2|f/|Ai,A3) 



(2S/vr2) L » pfP (g' « N) 

A^SeF, L<^pfP (g^^N) ^^^^ 
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C. Uncorrelated case 



Further calculations are different for GCE and CE. We start from the most simple GCE 
case. In this situation the averaged sum over the three different states Ai, A2, A3 in Eq. (pSf ) 
is essentially the integral over three energies, corresponding to these states, multiplied by 
the normalizing factor A"'^ and the four-point correlation function R4, that is responsible 
for the level repulsion for small energy differences: 

7(e) = 27rU^A-' J dtide2de^f2{l - /i)(l - h)6{t + - ei - e3)i?4(e, ei, £2, £3) (43) 

Now we set e = 0. As was already mentioned, an analysis gives the same energy dependence 
7(e, T = 0) as the temperature one 7(e = 0,T) obtained below, apart from the numerical 
coefficients. 

The characteristic scales of variation for the Fermi functions and the correlation function 
i?4 are T and A respectively. So for T ^ A the correlation function can be replaced with its 
asymptotic expressions for large values of arguments, i.e. i?4 = 1 [^. In physical terms this 
means that the level correlation does not play any role for the electron-electron scattering 
rate provided T S> A. Note also that in this limiting case the GCE and CE situations 
coincide essentially, for the number of excited quasiparticles is large. 

Direct calculation of the integral gives for T ^ A: 



luc- 2 ^3 



2 



(S/vr)(T2A/E,2) L » {g^ « N) 

(7rVl6)(TVe^) L<^pFf {g^^N) ^^^> 



(Fig. 4). The appearance of the mean level spacing A in the expression for the case where 
the level correlation is absent, should not be misleading: it just stands for a combination 
(z/L)~^, and is introduced for convenience. 

The upper line of Eq. (H3) (range III on Fig. 4) corresponds to the case of "large" dots. 



The dependence 7 oc is universal and is not sensitive to the geometry of the sample. 
However, the coefficient depends both of the sample's size and the mean free path. It is 
valid for T <^ Ec-, and gives in this parameter range 7 -C A. Hence the spectrum is discrete 
and this result can not be obtained in the perturbation theory. Note, however, that the if 
one formally applies the perturbative, e.g. those, developed in [|] or |TB|, methods to the 
range T <C Ec-, one obtains the same result 7 oc T^A/E"^, 7 ^ A. We have given above a 
rigorous, self- consistent derivation of this result. 

The lower line of Eq. (|^) is essentially the same result that appears in the clean limit 
for bulk 3D system. It does not contain either the size of the system or the diffusion 
properties (such as the mean free path). However, the range of validity for this result 
(L ^ pfP,A <^ T <^ Ec) is rather different. This range includes both cases discrete and 
continuous spectra, and consequently we have an overlap between the perturbation theory 
and the RMT calculation (region between curves 1 and 3, Fig. 4). Note that this "clean" 
behavior is observed in small enough dots even for small temperatures Tr <^ (pi?/)~^ (cf 

m 
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D. Effect of level correlation 



For T -C A the correlation function can be replaced by its expansion for small arguments; 
in the particular case of the Gaussian Unitary Ensemble (GUE) one has (see Appendix B): 

i?4(0, uj,uj + n,n) = (7r/A)^2(2i2625)'^(cJ®r]^ - 2cjW + u^n^) 

Straightforward calculation gives: 

W - (f ) 2., (f ) " - 1^ - 1000 (45) 



In particular, one obtains 

'Jcorr ~ 



/5S(ATV£;2)(7rT/A)^2^ L > {g^ <^ N) 
(7rV8)(TVe^)(7rr/A)i2, L « p^/^ (^2 » ^) 



(46) 



The twelfth power in the result can be easily explained. One needs to find four energies 
close to each other; these four energies form six pairs, and the contribution of each pair is 
proportional to (tu/A)^ in the GUE, u being the difference between the energies in a pair. 
Consequently one obtains the extra factor (T/A)^^ in comparison with the uncorrelated case. 
In a similar way, the contribution from each pair is proportional to |a;/A| in the Gaussian 
Orthogonal Ensemble (GOE) and to (cu/A)'* in the Gaussian Symplectic Ensemble (GSE), 
and so the results are 7 oc -fuc{T/A)^ (GOE) and 7 oc 7„c(T/A)24 (qSE). 

The CE case is more complicated. As we argue above, for T ^ A it yields the same 
results as the GCE, and so the only interesting situation is T <C A. In the expression ( pSf ) 
the functions / are still Fermi distribution functions, however the chemical potential needs 
now to be pinned to one of the energy levels: = ei + 0. This fact changes all results 
considerably, as the probability to find fi in some gap between levels depends no longer on 
the width of this gap p6| , ^ . The arbitrariness of selection of the pinned level is removed 
by the averaging with a weight function that is centered around the "mean value" of the 
chemical potential /i and has the support Sfi: A -C 5/i -C /i. If this weight function is chosen 
to be a step function, after taking a limit — one obtains: 

7A = 27rf/2A 5(eA + eA,-eA, -eA3)5(ez-/i + 0)/A,(l-/Aj(l-/A3) (47) 

In principle, the summation over / contains terms with / = Aj. These particular terms are 
reduced after the disorder averaging to the integrals over three energies with the four-point 
correlation function i?4. But upon setting e = these terms vanish just because of the 
correlation function, containing two equal energies (e; = 0). However for non-zero e these 
terms yield the results: 



7(e,T) 



7nc(e/T)4(r/A)i2 GUE 
7nc(e/T)2(T/A)6 GOE 



(48) 



All other terms contain the five-point correlation function and hence can be omitted. For- 
mally for e = one obtains, e.g., in GUE 7ce ~ 7uc(2^/A)^°. 
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V. CONCLUSIONS 



In conclusion, we have investigated the out-scattering time appearing due to the 
electron-electron interaction. Two different cases should be distinguished: continuous spec- 
trum (7 ^ A), where perturbation theory can be applied, and non-perturbative case of 
discrete spectrum (7 ^ A). 

In the case of continuous spectrum the out-scattering time is essentially the same as 
the inelastic scattering time entering the problems of universal conductance fluctuations 
and persistent currents. In 3D and OD systems it coincides also with the phase break- 
ing time, while in 2D and quasi-lD cases these times differ considerably. In 2D case in 
some range of parameters we have recovered the earlier results jl^^l^, but we interpret 
it as the out-scattering time rather than the phase-breaking time. Also an intermediate 
parameter range between 2D and OD systems is investigated. For quasi- ID systems we have 
obtained principally the results (|TBp-(^). In purely one-dimensional case for large enough 
temperatures Tr ^ {ppaY out-scattering time occurs to be proportional to T~'^/^ as 
well as r^, however the former occurs to be considerably shorter. For lower temperatures 
{pFCi)~'^ -C Tt <^ {pFa)~'^ the out-scattering rate is proportional to T^/"^, and becomes 
larger than the temperature. Also a transitional region between ID and OD system exists, 
and we have investigated different regimes of the diffusion. In particular, for T ^ Ec the 
spectrum is always continuous, and close enough to E^. we obtain zero-dimensional behav- 
ior: 7 ~ oc T. Ranges of parameters corresponding to different regimes are displayed on 
Fig. 3. 

For temperatures below the Thouless energy the spectrum is discrete. In this case the 
out-scattering rate coincides with the phase-breaking one, and it is reasonable to speak just 
of the inelastic scattering rate. We have shown that for large enough systems L ^ p^pl 
{N ~ ei^/A > the latter behaves itself as 7 ~ T'^A/E^ (A < T < EJ, while in small 
systems N <^ g"^ the universal dependence 7 ~ T"^ /ep is found. For T <^ A we have obtained 
an abundance of results for small - large / CE - GCE / GUE - GOE cases. In spite of 
this abundance one should clearly understand that the inelastic scattering rate due to the 
electron-electron collisions is in this parameter range vanishingly small. So for real systems 
it is necessary to look for another mechanisms of inelastic scattering. Electron-phonon 
scattering seems to be a good candidate; it has been recently discussed for mesoscopic 



systems |^8[- Also a coupling to the environment produces an inelastic rate in the nearly- 
closed system. 
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C. M. Marcus, A. D. Mirlin, H. Schoeller, and U. Sivan are also gratefully acknowledged. 
The work was supported by the Alexander von Humboldt Foundation and SFB195 der 
Deutschen Forschungsgemeinschaft. 

Note added. As this paper was being prepared, a number of related works came to my 
attention. First, model similar to the "discrete spectrum" situation in quantum dots have 
been recently studied with the help of the random band matrices, and the results are rather 
similar Then, Kamenev and Gefen studied the role of the external enviroment 



in the inelastic broadening, and found this effect to be very strong. They relate this fact 
to the phenomenon of the Coulomb blockade. Finally, the recent unpublished results by 
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Altshuler, Gefen, Kamenev and Levitov |^ (AGKL) and the comment on their work by 
Imry |32] are so close to the results obtained above that some discussion is required. AGKL 
study the same problem, and reproduce, in fact, the matrix element (^Tl). They interpret 
it as an overlap between one-particle |A) and three-particle |Ai, A2, A3) states. Consequently 
this overlap should be compared with the three-particle level spacing A3 ~ A^/e^, and this 
comparison creates a new energy scale E* = (EcAy^"^. For E > E* the three-particle states 
are well mixed by the Coulomb interaction, and the broadened peaks (which are resolved for 
E < Ec) are essentially a mixture of many-particle states. On the other hand, for E < E* 
the single-particle state is mixed with one three-particle state, one five-particle state, and so 
on. In this sense the inelastic scattering rate is zero: the state does not decay at all. AGKL 
describe this situation as an analog of the localization transition on the Bethe lattice. I do 
not want to address this problem here; however, the results obtained above for E < E* can 
be interpreted rather as the width of the "envelope", formed by the many-particle states 
around the single-particle one. I am indebted to the authors of all these papers for the 
possibility to become acquainted with their results prior to publication. 



APPENDIX A. DEPHASING AND QID - OD CROSSOVER 

Below results for the phase-breaking time in the range of parameters intermediate be- 
tween quasi-one- and zero-dimensional case are obtained. In principle, one has to perform 
calculations similar to those of Ref. [|1^] in restricted geometry, and this does not look quite 
hopeful. However, as we are interesting in the result up to the numerical factor only, it is 



reasonable to use the method developed by Stern, Aharonov, and Imry |o^;0^ that was later 
applied to a calculation of the quasiparticle lifetime in a quantum dot ||18| . In this approach 
the phase uncertainty -P(to) accumulating by the electron due to the interaction with the 
environment (in our particular case it interacts with the electromagnetic fluctuations) is 
calculated; the time to when this phase uncertainty becomes of order unity (we will set it 
exactly unity) is associated with the electron lifetime (phase-breaking time). This phase 
uncertainty is given by [|18 



^(^0) = / dtdt' / du coth — ^ — Im e*-(*-* ) {S, + 52 + ^3 + ^4) 

(Al) 



with 



5*1 = exp{iq[xi{t) - Xi(t')]), S2 = exp(iq[x2(t) - X2it')]), 

S3 = exp(iq[xi(t) - X2(t')]), 'S'4 = exp(iq[a;2(t) - Xi{t')]) {A2) 

In these expressions (to be the direct generalization of those of Ref. for the case of finite 
temperatures) xi{t) and X2{t) are arbitrary (quasiclassical) electron paths, and the averaging 
over these paths is supposed; e{q^,uj) is the dielectric susceptibility. The expression (Al) is 
derived for an infinite system, however, as it varies on scales of order of the elastic mean 
free path / ^ a, this expression is valid in the diffusive regime, and the bulk results may 
be substituted for the dielectric susceptibility. It is also important that the q = mode 
does not contribute to the sum (Al): it is absent in the open systems while in the closed 
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ones the corresponding contribution is forbidden by the charge neutrahty. So the difference 
between open and closed systems is not quite important, and the system is assumed to be 
closed: q = 7r(na./L, n^/a, n^/a), rix + ny + nz > 0. The results for the open systems depend 
on boundary conditions, but in all cases differ by a numerical factor of order unity only. 

After the averaging over paths is carried out in exactly the same way as in Ref . [|18| , one 
obtains: 



48T D 



E 



1 



/ rft+ / — sin(2Tt)x 
Jo Jo t 



^^nl + {L/annl + nl 



X exp 



-TT 



'■E,t {nl + [L/a] 



2/ 2 I 2 

{riy + 



Here D is the three-dimensional diffusion coefficient and we have taken into account that 
the main contribution to the integral over frequencies comes from the region \uj\ < 2T. The 
summation is restricted by the condition q <^ I. For D /a? ^ T (this condition excludes 3D 
situation) only the values = = are important in the summation (A3). 

case and the result is given by Eq. 



In the hmiting case Ec ^ (^^e inverse situation corresponds to the "true" quasi- ID 



the integral can be easily calculated. One gets: 



Pito) 



AST D 



TipplLa? E, 



c n=l 



to arctan 



2T 



2T 



(A4) 



The phase-breaking time is defined as the time to when the phase uncertainty P{to) 
is equal to unity. 

The case Ec^ T corresponds to the "true" OD situation. One has 



96T 



n^ppLa'^ El 



I 



(A5) 



This is the result by Sivan, Imry and Aronov [T^ up to the numerical factor. It is seen, 
however, that in this region the spectrum occurs to be discrete: r^^ ^ A. Thus, the 
OD region is not subject to our analysis and should be treated by another methods. (See, 
however, the discussion after Eq. (PD). 

The case r^^ -C E'c <^ T is intermediate between quasi-lD and OD regimes. Only terms 
with E^n'^ < T are important in the summation (A4), however due to the condition T ^ Ec 
the summation can be extended to the infinity. One obtains: 



1 



that is a new result. In the whole region <^ Ec <^ T the spectrum turns out to be 
continuous and the result (A6) is valid. 
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APPENDIX B. FOUR-POINT CORRELATION FUNCTION IN CUE 



Below we follow the generalities given in Refs. |]2D|. As an explicit expression for the 
four-point correlation is not given anywhere in literature to the best of our knowledge, we 
derive it for the most simple GUE case. It is convenient to use the dimensionless energies 
X = vre/A throughout this Appendix. 

The first step is to define the functions 

Yi{xi, X2, . . . ,Xi) = Y^ S12S23 ■■■Sii (Bl) 
p 

Here 

sinx 



s 



s{x) = , {B2) 



X 



and the summation is carried out over {i — 1)\ different permutations of indices. Hence, 

Yi{x) = l; Y2{xi,X2) = sl^; Y3{xi,X2,X3) = 2S12S23S31; 

3^1? ^2, X3, X4) — 2S12S23S34S41 + 2S13S34S42S21 + 2S14S42S23S31 

Now the correlation functions be expressed as 

Ri = 1; i?2(a;i,X2) = -Y2{xi,X2) + Ri{xi)R2{x2); 

R3{xi,X2,X3) = Y3(xi,X2,X3) + Ri{xi)R2{x2, X3) + Ri{x2)R2ixi, X3) + 

+Ri{x3)R2{xi, X2) - 2Ri{xi)Ri{x2)Ri{x3); 



Ri{Xi,X2, X3, X4) = -Yi{xi, X2, X3, X4) + {i?i(Xi)i?3(x2, X3, X4) + Ri{x2)R3iXi,X3, X4) + 
+Rl{x3)R3{Xi,X2, X4)+Rl{x4)R3{xi, X2, 0:3) } + {i?2 , X2)R2{X3, X4)+R2{XI,X3)R2{X2, X4) + 
+R2{Xi, X4)r2{x2, X3)} - 2{R2{Xi, X2)Rl{x3)Ri{x4) + -R2 (^1 , X3)i?i (X2)-Rl (X4) + 
+R2{Xi, Xa)Ri{x2)Ri{x3) + R2{X2-, X3) Ri{xi) Ri{Xi) + R2{X2-, Xa)Ri{Xi)Ri{x3) + 
+R2{X3, X4)Rl{xi)Ri{x2)} + 6Ri{xi)Ri{x2)Rl{x3)Ri{x4) 

After some algebra one obtains an explicit expression for the correlation function R4: 

R4{Xl, X2, X3, X4) = 1 - 2{S12S23S34S41 + S13S34S42S2I + S14S42S23S3I } + 
{^12^34: + "5l3"S24 + "5l4'^23} + 2{Si2S23S3i + S12S24S41 + S13S34S41 + S23S34S42} — 

-{S?2 + 43 + + ^23 + ^24 + S34} (^3) 

In our particular case Xi + X3 = X2 + x^, and, taking into account that the correla- 
tion function depends on three differences of arguments only, one arrives to an expression 
i?4(0, x,x + y,y). From (B3) after cumbersome calculations one obtains an expansion of R4 
for x,y <ti 1: 

i?4(0, x,x + y,y) = - 2x%^ + xV) (^4) 
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The 12th power can be easily explained. Pair correlation function for small arguments is 
proportional in GUE R2{x) oc x^. As one has in six pairs of arguments close to each 
other and each pair produces the second power, the total expansion starts from the power 
twelve. In GOE R2{x) oc x, and we may easily conclude that the analogous expansion starts 
from the sixth power. 

As could be expected, for a large value of the arguments x, y ^ 1 all functions Sij 
are small, and R4 = 1. This fact means just that the levels are uncorrelated on distances 
exceeding the mean level spacing A. 
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FIGURE CAPTIONS 

Fig. 1. The diagram representation of the diffusion propagator. 

Fig. 2. The diagrams contributing to the interaction block (crossed rectangle on Fig. 1) 
for the UCF diffusion propagator. Wavy hnes and dashed hnes with a cross denote the 
electron-electron interaction and bare impurity scattering (tti^t)"^ respectively. 

Fig. 3. Different regimes for fD - OD crossover. Curves: 1: g = Ec/A = 1; 2: Tr — 
{l/Lfippay (or 7 = E^ippay); 3: Tr = {l/Lf{ppay (or E, = T{pFa)-^); 4: Tr = 
{1/ LYippaY (or EcT^ = 1); 5: Tr = {1/ L)^ (or T = Ec- Regimes I - V are described in Sec. 
IDD. 

Fig. 4. 3D - OD crossover. Curves: 1: T = E^; 2: T = A. Line 3 (7 = A) separates 
the cases of continuous and discrete spectrum for L ^ p^t^. Regimes: I - clean limit, 
7 ~ T^/ci?; 11 - bulk diffusive limit, 7 oc T"^/^; 111 - zero-dimensional diffusive limit; in the 
ranges IV and V the level correlation is important. 
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